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2 nd Order Parallel Splitting Methods for Heat Equation 
Aziz, M. a and Rehman, M. A. b 

Department of Mathematics 

University of Management and Technology, Lahore, Pakistan. 
azizkb9700@gmail.conT, aziz3037@yahoo.com b 

Abstract 

In this paper, heat equation in two dimensions with non local boundary condition 
is solved numerically by 2 nd order parallel splitting technique. This technique 
used to approximate spatial derivative and a matrix exponential function is 
replaced by a rational approximation. Simpson’s 1/3 rule is also used to 
approximate the non local boundary condition. The results of numerical 
experiments are checked and compared with the exact solution, as well as with 
the results already existed in the literature and found to be highly accurate. 

Keywords: heat equation, 2 nd order numerical methods, method of lines, parallel 
algorithm 

Introduction 

Several authors have been studied the initial- boundary-value problems (pdes). 
These equations are mathematical models for real world phenomena which often 
take the form of equation relating to various quantities and their derivatives. The 
motion of particle in the straight line, motion of a missile (projectile motion), 
function of a nuclear reactor, transfer of heat, vibration of a particle and chemical 
reaction etc. All such problems are governed by elliptic, hyperbolic or parabolic 
partial differential equations which may or may not be homogenous, in one, two 
or three dimensions with non local boundary conditions, along with initial 
conditions existing in the literature. In real life problems, parabolic equations with 
an integral boundary condition(s) have wide applications, and some time we are 
required only their numerical solutions. The concentration of a chemical reaction 
conveniently by employing absorption of light at suitable frequencies. The 
summation property of electronic and limited availability of light beam, the 
computational process of used chemical gives rise to quadrature of concentration 
in spatial region from which the light beam pass through is described by a 
parabolic equation. We can also use this method for smoothening of data and to 
calculate time integral of spatial term computation data relating to kinetic rates of 
chemical reaction. Thus parabolic partial differential equations with nonlocal 
boundary condition(s) have great impact in diverse fields such as 
electrochemistry, heat transition, biological and medical sciences and population 
dynamics. Mostly authors used explicit or Crank-Nikolson finite difference 
schemes for the numerical solution of these partial differential equations. In 1994 
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B. J. Noye. et al. have worked out a numerical scheme for the solution of heat 
equation along with initial condition, boundary conditions and non-local boundary 
specifications and they used sufficiently smooth data which satisfy the required 

compatibility conditions in the scheme. They used two different explicit schemes 

1 

of finite difference which is of second order but for S = - this scheme is of forth 

6 

order. They have employed the well-known two finite difference scheme 1 st is [1], 
[5] FTCS and 2 nd is [1], [9] FTCS, and the integral was approximated by 
Simpson’s 1/3 rule to solve the heat equation. He reduced the time consumed 
4330.7 sec to 410.3 sec as well as the results were much better than the previous 
results existing in the literature [2]. 

The same problem was also discussed by A.B. Gumel. et al.. They 
replaced the space derivative by central difference approximants of 2 nd order and 
approximated the quadrature existing in the recurrence relation by the Trapezoidal 
rule [3] and the matrix exponential in the recurrence relation by Pade’s 
approximant [4]. They suggested a different scheme to solve the heat equation. [3] 
used the method of line semi discretization and transformed model PDE into the 
system of linear, 1 st order ODE’s, and this system of ODE’s satisfied a certain 
recurrence relation. In solution of ODE a matrix exponential term was produced 
which was approximated with pade approximation, and for unknown co(t), he 
used Trapezoidal rule. [3] also presented the result in relative error form and 
compared with [1], [5] FTCS explicit method . Yunkang Liu worked on heat 
equation in one dimension with non local boundary conditions. Several authors 
have discussed PDE’s in one dimension as well as in two dimensions with non 
local boundary conditions [5]. In 1996 Twizell and Taj introduced a new rational 
approximation to the matrix exponential function, which creates while solving 
parabolic partial differential equations. In 1997 Taj and Twizell extends their 
work to fourth-order parallel splitting methods having all properties as in the third 
order but with much batter error of approximation as compared to previous 
results. These methods do not require any complex arithmetics and are third-order 
accurate in space and time as well as L-acceptable, tested and found more 
accurate than the results existing in the literature [9]-[ 10]. In 2003 Mehdi 
Dehghan started working on the theory of parallel splitting method and applied 
this technique for the solution of hyperbolic and parabolic PDEs. Most of his 
effort appeared from 2003 to 2005 and he used different methods to get a hold for 
approximate solutions of the heat equation [8]. In 2008 Jurgen Geiser applied 
fourth order parallel splitting method for the solution of heat equation by using 
ADI (Alternating Direction Implicit) methods and LOD (Locally One- 
Dimensional) methods, which are standard splitting methods of lower order, e.g. 
second order. He extended this second order to fourth order based on locally 
higher-order LOD methods. In 2007 the heat equation with boundary integral 
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specification and mixed order parallel splitting method was discussed by M. A. 
Rehman. For the numerical solution of non classical PDEs, S. Ali Mardan and M. 
A. Rehman have introduced a hydride method. 

We shall use Simpson’s 1/3 rule for approximating the integral condition 
i.e. non-local boundary condition. One way of solving parabolic partial 
differential equations is by the use of method of lines (MOL), which transform 
parabolic partial differential equations into a system of ordinary differential 
equations which can be written in matrix-vector form as 

® = 51/(0+ <K0 (1) 

Where 5 is a square matrix, <p(t) results from non-homogeneous 
boundary conditions and U(t) is the solution vector at time t. 

Development of the Method 


Consider two-dimensional heat equation 

du(x,y,t ) d 2 u(x,y,t ) 

dt dx 2 

0 < x,y < 1, t > 0 
Subject to the given initial condition 


+ 


d 2 u(x, y, t) 
dy 2 
( 2 ) 


u(x,y, 0) = g(x,y ), 
and the boundary conditions 

0 < x, y < 1 

(3) 

u(0,y,t) = f 0 (y,t), 

t > 0 

(4) 

u(l,y,t) = A (y,t) , 

t > 0 

(5) 

u(x, l,t) = g^x.t) , 

t > 0 

(6) 

u(x, 0,0 = floWffl(t) 

, t > 0 

(7) 


With nonlocal boundary condition 


n 

II 

oo 


u(x,y, t)dxdy — g{t) 0 < x,y < 1 (8) 


Where 

fo (y> 0, fi (y, 0, g 0 O- 0, 9i (*, 0 and 

fu(t) are functions which are continuous and known, where as u and co are 
the unknown functions to be determined. 


Discretization 


The intervals 0 <x < 1 and 0 < y < 1 are divided into N + 1 subintervals of 
1 

equal length h and h — and t is the time variable and length of each time step 

is equal and is /, and occupying space R = [0 < x < 1] x [0<y< 1] x [t > 0J 
gives a cuboidal mesh of cuboids of equal volume h 2 xl and the co-ordinates of 
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vertices of each cuboids are ( x m , y m , t n ) = (mh , mh , nl ) where m = 
0,1,2,3, ...,1V +1 an d n = 0,1,2,3,... The boundary dR of whole space is 
consisting of linesx = 0,x = l,y — 0,y — 1 and t — 0. 


Derivation 

The space derivative 9 and 


d 2 u(x,y,t ) 
dy 2 


are approximated as under 


9 -^^- = ± {u( x - h,y, t) - 2 u ( x , t) + u(x + h,y,t)} + 0(h 3 ) (9) 

and 

9 {u(x,y — h,t) — 2u(x,t) + u(x,y + h,t)} + 

0(h 3 ) (10) 

Applying these approximations to the each mesh points of the grid, at time 
level t = t n , we get a system of N 2 first-order ordinary differential equations. This 
may be written in matrix form as 


dU(t) 

—^- = BU(t)+<p(t), t> 0 

With initial distribution 

1/(0) = g(x,y) 

The matrix B is the sum of two square matrices B\ and Ih of order N 2 . 


5 .=tt 


h~A 


h A 


O 


O 


h 2 A 

h 2 A 


where A is given by 


-2 1 
1 -2 1 

1 

O 


O 


-2 1 
1 -2 
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and 


2 h 2 


-21 I 
I -21 I 


O 


I -21 I 

O 1 -21 

Where both B i and B 2 are of order N 2 x N 2 and I is the identity matrix of order 
N. The vector (pit), sum of two vectors (pi(t) and (piil) of order N 2 x 1 arises from 
the use of the boundary conditions in above two approximation. 

Then solution U (t) satisfies the recurrence relation 

U(t + l) — exp (IB) U{t ) + 

t+i 

J exp ((t + l — s)B)cp(s) ds; t = 0,1,21, ... 

t 

Where / is constant time step in the descretization. To approximate the matrix 
exponential function in a new rational approximate, for a real scalar 9, given by 


E 2 {6 ) - 


1 + (1 — cl)9 + 

1 — aO + (a — i)0 2 


The unknown U n will be ordered in the form 

In which (t) = [l/ lfl (t), U 21 (t), U 31 (t),... , U N1 {t), ... , U 1N (t), U 2 N it),... , 

9 = V9i,i’92,i>93,i> ■■■ <9 n,i< ■■■ < 9i,n >92,n> ■■■ >9n,n] > 
and the solution of system satisfy the recurrence relation (2). Where 

t+i 

j exp ((t + l-s)B) v ( s )ds=W lV (s 1 ) + W M s 2 ) 


Where =£ s 2 and weights VKj(t = 1,2,) are matrices. It can easily be shown 
that when 

(pis) = [1,1,1,. . . ,1] T , then 


where 
and when 


Wi+W 2 =M 1 

M, ={B l (exp(lB)-I)} 

<p(s) =[s,s,s,. . . ,s] T , then 
Si Wi+S 2 W 2 =M 2 
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where M, = B' {(exp(/B) ■(/ I III I ft (exp (IB) - 1 )[ 

Using 0=(S in E 2 (fl)= 

exp(ZS) = PQ we have 
P(Zfl) = (/ - alB + ((a - ^Z 2 # 2 )- 1 
and 

Q(Zfi) = 7 + (1 - a) IB 

Taking s x = t , s 2 = t + Z 

and then solving and replacing by exp(ZB) gives 

U(t + l) = exp(Z5)Z7 (t) + 

l -[G{lB)(p{t) + K(lB)<p(t + /)]; t = 0, l, 21,... 


In which, G(IB) = I - alB + (a - —)(//?)" 


1 2 1 

K(IB) = I - alB + (a - -)(ZP) 2 I-2(a--)IB 

Parallel Algorithm 

In parallel Algorithm exp(Zfi ), G(IB) and ( IB ) are decomposed in their partial 
fraction form given below 


Where 


exp(Zfi) = ^ q t (7 - AJB) 1 

/=! 

Where 

1-a + A 1-a + A, 

q, =- L , q 2 =-- 

1 2 


G{lB) = Y j q i+2 {I-^lB) 1 


_ ji _ 

Km=i,q M v-4iBr‘ 
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where 


Ay — A, 

Where A : and A 2 are given below 


1 - 2a + A 1 - 2a + A, 

4 5 =—;- :—^6 = 


A, ~ Ay 


and 


Hence equation 


A = 


2a -1 


4 = 


a + 4a 2 - 4a + 2 
2a-1 


a 


- \[a — 4a + 2 


U(t + l) = exp(lB)U(t) + 

l -[G(lB)cp(t) + K(lB)q>(t + l)\,t = 0, /, 21,... 


become 


U(t + l) = YAT\ qi U R{t)+ l - ( q i+1 (p{t ) + q i+4 (p(t +1 )) 
»—1 ^ 


Where 

Hence 


A. = (I-AlBy\i = 1,2 


U(t + l ) = 

i=l 


Where y* are the solutions of the system. 

/ 


A y,=<i,u { t ) + ^[ ( ii + 2 ( p{ t ) +( ii + A ( p{ t+l )} 


Numerical Experiment 


u(x,y,t) = e x+y 
u(0,y,t) =e y+2t 
u(l,y,t) - e l+y+2t 
u(x 9 l,t) = e 
u(x,0,t) = g o (x)a>(t) 
g 0 (x) = e x 
fi(t) - (e-\) 2 e 21 
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Oj(t) = e 1 ' 


The analytical solution of the problem is 

u(x, y,t) = e 


x+y+2t 


Table 1. 

1 

Relative errors for n = 49 and Z =-when co (!) is given 

15000 w b 


Node points 

x = y 

Theoretical 

solution 

2nd order 

finite difference [3] 

New Scheme 

0:1 

9.0250 

2:5562 x 10' 7 

1.171 x 10- 6 

0:2 

11.0232 

2:8302 x 10' 6 

2.881 x 10- 6 

0:3 

13.4637 

4:1114 x lO' 6 

4.151 x 10' 6 

0:4 

16.4446 

4:6798 x 10' 6 

4.714 x 10' 6 

0:5 

20.0855 

4:5541 x 10' 6 

4.586 x lO' 6 

0:6 

24.5325 

3:8804 x 10' 6 

3.913 x 10' 6 

0:7 

29.9641 

2:8540 x 10' 6 

2.890 x 10' 6 

0:8 

36.5982 

1:6804 x lO' 6 

1.722 x 10' 6 

0:9 

44.7012 

1:0154 x lO' 6 

6.340 x 10' 7 


The Treatment of Non Local Boundary Condition 

By using Simpson’s 1/3 Rule. 

II j2 M M 

11 u(x, y, t)dxdy * — 

0 0 ^ i =o ;=0 

Where M — N + 1 

^0 = = 1 

Cj = 4 for odd i and c t = 2 for even Z 
Above equation may be written as 

r r h 2 ^ 

J J u(x, y, t)dxdy » — 2^c,U ifl + ^(0 

0 0 ^ 1=0 


Where 


h 


2MM 




i=0 7=1 


Integrating the boundary condition with respect to x it gives 

i i 

Jm(x,0,0 = co(t)j g 0 (x)dx 


Once again use Simpsons’ Rule we get 
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4 i=0 n 


<w(t) = 


m( 0 - /?(t) 


fa 1 /:* ( x ) dx ) 


Table2. 

1 

Relative errors for n = 49 and Z =- When <u (/) is unknown 

_ 15000 ___ 


Node points 

v = y 

Theoretical 

solution 

2nd order finite 
difference [3] 

New Scheme with 
Simpson’s rule 

0:1 

9.0250 

3/3994 x 10 4 

6.049x 10 6 

0/2 

11.0232 

3/0427 x 10 ' 6 

3.207 x 10 ' 6 

0/3 

13.4637 

2/6002 x 10 4 

5.420x 10 ' 7 

0/4 

16.4446 

1/8408 x 10 4 

1.422x 10 ~ 6 

0/5 

20.0855 

1/1595 x 10 4 

2.512x 10- 6 

0/6 

24.5325 

6/3782 x 10 5 

2.764x 10' 6 

0/7 

29.9641 

2/9338 x 10 5 

2.349x 10 ' 6 

0/8 

36.5982 

5/1982 x 10 6 

1.526x 10 6 

0/9 

44.7012 

3/3960 x 10 6 

5.940 x 10 ' 7 


Conclusion 


Using given algorithm the problems are solved for n = 49 and l — igooQ . 

We have developed 2 nd order parallel splitting method for heat equation in two 
dimensions and non local boundary condition was approximated by Simpson’s 
1/3 rule and exponential term in recurring relation by new rational approximate. 
First we solved the problem when the value of co (t) was given and then for 
unknown co (t). The method is Lo stable and used for two dimensional heat 
equation with non local boundary conditions. The results for relative error 
obtained by using above numerical method compared with the results [3], which 
are better. It is also observed that there is no oscillation over the entire interval 
and the method behaved smoothly. 
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Generalization of TOPSIS from Soft Set to Fuzzy Soft Sets 
in Decision Making Problem 

Saeed, M., Anam, Z., Kanwal,T. Saba, I. Memoona, F. and Tabassum, M. F.* 

Department of Mathematics, School of Sciences, University of Management and 

Technology, Lahore, Pakistan 

* farhanuetl2@gmail.com 

Abstract 

As of late, the issue of basic decision making has discovered imperative centrality. 
It has procured centred significance particularly for the issues identified with 
incorrect environment. The technique for order of preference by similarities to 
ideal solution (TOPSIS) is the multi criteria choice examining technique utilized 
for the choice. We apply the generalized result of TOPSIS on soft set to TOPSIS 
on fuzzy-soft-set. 

Keywords: TOPSIS, Soft-sets (SS), fuzzy-soft-sets (FSS), weighted-normalized- 
decision-matrix (WNDM) andfuzzy-sets (FS). 

Introduction 

In real life, we face many problems of decision making in different domains of 
sciences such as engineering, medical science and economics. To address 
uncertain situation and treat ambiguous data, Zadeh [1] presented the idea of 
fuzzy sets in 1965. Introduction of fuzzy set theory revolutionized the entire 
mathematical sciences. Many researchers contributed in the development of fuzzy 
sets and its applications in various fields of life. 

In 1980 Thomas Saaty [2] introduced Analytical Hierarchy process. These 
are most useful techniques for dealing difficult decision making and help the 
decision makers to make accurate decisions. This theory consists of many 
applications such as, medicine, computer science, control engineering and 
artificial intelligence etc. 

To ease this situation of decision making a methodology was developed 
by Hwang and Yoon [3] in 1981, named as TOPSIS. But TOPSIS on soft set is 
useful for discrete situations. In 1999 Molodtsov proposed the theory of soft-sets 
as a generalization of fuzzy set, which can efficiently handle a number of 
parameters simultaneously. 

The idea of fuzzy-soft-set was first presented by Maji et al. in 2001. In 
2002, firstly the technique of fuzzy-soft-set was used in decision making. In 2013, 
Maji et al. presented the idea of the neutrosophic-soft-set. For decision making 
many techniques are used such as soft expert set, AHP, interval valued fuzzy soft 
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matrix, TOPSIS etc. These techniques are proved to be very helpful in decision 
making. M. Saeed, Sana. A. and N. Rubi [4] compared the two different technique 
, fuzzy-soft-expert-set and AHP. The computational cost of Fuzzy-Soft-Expert- 
Set system is minimum than AHP procedure for the same outcome. They 
conclude that AHP and fuzzy-soft-expert-set gives the same result. 

Zulqarnain. M and M. Saeed [5] have applied fuzzy-soft-set in decision 
making. The idea of fuzzy-soft-set has been presented earlier [6]. But decision 
making on fuzzy-soft-set by applying TOPSIS is presented for the first time in 
this paper. In this paper, we generalized the concept of TOPSIS from soft-set to 
fuzzy-soft-sets in decision making. In this paper the author present some basic 
concepts associated to soft-set and fuzzy-soft-set and procedure of TOPSIS. Also 
develop the methodology of decision making via TOPSIS on fuzzy-soft-set and 
conclude the results. 

2. Preliminaries 

In this subcategory, we present the fundamental concepts and outcomes of 
soft-set-theory, which would be beneficial for additional dialogues. Maximum 
descriptions and outcomes obtainable in this subdivision may be established in [7- 

9]. 

Definitionl 

“Let S be a set of parameters and IT be an initial universe set. Let IT be the 
power set denoted by P (U) the and AcS. A pair (G, A) is called a soft set over IT, 
given by 

G: A —> P (Uf 

Where G is a mapping [10] 

Fuzzy-Soft-Set in Decision-Making 

In this section, we will discuss some elementary description of fuzzy soft 
set (LSS) and some outcomes, which we use in further debate. Most of them are 
originating in [ 11,12]. 

The parameter space S may be written as 

S 2 [5, .ufi }. 

Let each parameter set Yj characterize the j th class of parameters and the 
elements of Bj represents a particular property set. Here we suppose that these 
property sets may be observed as fuzzy-sets (LS). 

Definition! 

“Let P (V) represents the set of complete fuzzy-set of V. Let XcL. A pair 
is ( Hi,Xi ) is called a fuzzy-soft-set over V. We may now describe a LSS ( Gj,Bj. ) 
which illustrates a set of objects having the parameter set [13], where Hi is a 
mapping given by H j : A, -> P(V) . [14] 
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TOPSIS 

The basic steps of TOPSIS are described in [15] 

Stagel. First determine the normalized-decision-matrix. 

The normalized value is denoted by x (/ is calculated as follows: 

x a = y<J|>S i =1 > 2 --« n and j = 1, 2,m. 

Stage 2. Evaluate the weighted normalized decision-matrix. The weighted 
normalized value Uij is calculated as follows: 

Uij = XijX Wi, 

i =1, 2 ..., n and j = 1, 2 ... m. [16] 


where Wi is the weight of the i th alternative and XT = \. 

i'=l 

Stage3. Calculate the best ideal solution ( B*) and worst ideal (B ) solutions [17, 
18], 

B = ((«“ Uij 1 • G D> (™ n Un 1 • e D ))=(w j»'=O,-,«} 

5" ~ h ni ' n m,, i«' e DM^Uu l,e D ( j) = 1 »=t 2 --") 

Stage4. By using the n-dimensional Euclidean space compute the distance. The 
parting procedures of each substitute from the best-ideal-solution and the 
worst-ideal-solution, respectively, as given below: [19, 20] 

I) //'>'•' i - 2 ." 


Stage5. Compute the relative closeness to the best-ideal-solution. The relative 
closeness of the alternative with respect X,- to B is defined as follows [21]: 


RC, = 

Stage 6. Rank the preference order. 


B, 


B+B■ 


,i = 1,2,.. 


n 


A Decision Making Method on Fuzzy Soft Set 

In this section, we deliberate the decision-making method by using TOPSIS 
on soft-set-theory. The detailed method, of each step, is shown below. 

Stepl. Defining the problem. Let us assume that 

DM = { D p ,p E I n } are the sets of decision makers, 
v — [vi i E I m ]denote the set of alternatives 
and x — [xj j E I n } is set of all parameters. Then FSS over v is a function defined 
by 


Hj : X t -> P(V) 
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where Hi a mapping is given by [22] 

Step2. Construct decision matrix D for each decision makers. 

Wher eD — V i j EIn d i j,d i j—fx i {Xj') is the criterion values of i th alternatives 
received from the criterion, X t is the parameter sets of decision makers D p and 
f x .is the soft set which was constructed by D p . 


v i 



*1 

d ll 

d 21 


v 


m 


d„ 


X 2 


d 12 
d 


22 


d 

d 


In 
2 n 


d„ 


d„ 



mxn 


Step3: Obtaining the WNDM of V. The WNDM is calculated as: 


y Imxn 

where v kt = Vi,/c,tG/ n 

Step4. 

n 

K u i) = ^ v ji 

i =1 

where, k (uf) is decision values of Uy.Thus the decision matrix of each alternative 
values for the deciders is expressed as 

R = {/cOi),. k(u n )} 

Step5. Ranking the preference order. 

An Application 

Stepl. Defining the problem 

Suppose that a car dealer has a set of various types of cars (universal set- 
alternatives) 

v = {v 1 ,v 2 ,v 3 ] 

which may be categorized by a set of all parameters 

Y = {y l ,y 2 ,y 3 }Forj= 1 , 2 , 3 . 

The parameters Y j , stand for "luxuries”, "automatic” and "manual” respectively. 


7 

11 

<N 

•• V m 

t 

21 

V 22 

'• v 2n 

ml 

V „,2 ' 

• • V 

mn 
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Assume that three decision-makers come to the car dealer to buy a car. 

Firstly, each decision-maker has their own choice about car and they 
consider their own set of parameters. 

Then they can construct their fuzzy-soft-sets. Next, by using the fuzzy-soft- 
set and TOPSIS-decision making method we select a car on the basis of 
parameters of decision makers. 

Step2. Construct decision matrix D for each decision- makers. We can 
construct fuzzy soft sets of decision-makers, in a tabular form respectively as 
Fuzzy soft sets of decision-maker, D 1 is 



+1 

33 

y 3 

V 1 

"0.3 

0.5 

0.2 

v 2 

0.2 

1 

0 

V 3 

0.4 

0 

0.3 

Fuzzy soft sets of decision-maker, 

D 2 is 




A 

? 2 

y 3 

v i 

"0.2 

0.5 

0.4 


0.3 

0 

2 

V 3 

0.1 

0.5 

0.4 

Fuzzy soft sets of decision-maker, 

D 3 is 




A 

? 2 

y 3 

v i 

"0.4 

0 

0.3 

v 2 

0.5 

0 

0.1 

V 3 

0.4 

0.2 

0.1 


Step3. Creating the weighted-normalized-decision-matrix V. 

Now compute the weights corresponding to each parameter. 

3 

v n= I fyxyi) /yi(yi)(vi)+/y 2 (yi)(vi)+/y 3 (yi)(v!) 

i=l 

= 0.3 + 0.2 + 0.4 = 0.9 

3 

V i2=^/yj( yl ) (. v 2 ) = /yi(yi)(v 2 )+/y 2 (yi)(v 2 )+/y 3 (yi)(v 2 ) 
i=i 

= 0.2 + 0.3 + 0.5 = 01 

3 

V i3=^/y i(yl) = /yi(yi)(v 3 )+/y 2 (yi)(v 3 )+ /y 3 (yi)(v 3 ) 

i =l 
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= 0.4+ 0.1 + 0.4 = 0.9 

3 

v 2l= ^/y i(y2) ( v l) = /yi(y2)(v 1 )+/y 2 (y2)(v 1 )+ /y 3 (y2)Oi) 
i=1 

= 0.5 + 0.5 + 0 = 01 

3 

v 22= /y i(y2) ^ 2 ) = fy 1 (y2)(y 2 )+fy 2 (y2)(y 2 )+ fy 3 (y2)(v 2 ) 

i=l 

= 0.0 + 0.0 + 01 = 01 

3 

v 23= ^ fyi{y 2 ) = fyi(y2)(v 3 )+fy 2 (y2)(y 3 )+ /y 3 (y2)(v 3 ) 

i =1 

= 0.0 + 0.5 + 0.2 = 0.7 

3 

v 3i= ^/y i(y3) ( v i) = /yi(y3)(v 1 )+/y 2 (y3)(v 1 )+ fy 3 (yi)(v{) 
i=l 

= 0.2 + 0.4 + 0.3 = 0.9 

3 

U 32 = I fyi(y 3 ) ^ 2 ) fyi(y3)(.v 2 )+fy 2 (y3)(y 2 )+fy 3 (y3)(y 2 ) 

i=l 

= 0.0 + 0.2 + 0.1 = 0.3 

3 

u 33= ^/y i(y3) (+ 3 ) = /y 1 (y3)(v 3 )+/y 2 (y3)(v 3 ) + /y 3 (y3)(v 3 ) 
i=l 

= 0.3 + 0.4+ 0.1 = 0.8 
Then the weight matrix is obtained as 

“0.9 1 0.9 

V= 0.1 1 0.7 

0.9 0.3 0.8 

Step4: Creating the decision matrix (vector), R. 

Now, calculate the individual elements of the R matrix. 

3 

k(y 1 ) = v u — v 11 + v 12 + v 13 — 0.9 + 0.1 + 0.9 = 2.8 

i =1 
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3 

k(v 2 ) = ^ V lt = v 21 + v 22 + v 2 3 = 01 + 01 + 0.3 = 2.3 

i=l 

3 

k(v 3 ) = ^ v lt = v 31 + v 32 + v 33 = 0.9 + 0.7 + 0.8 = 2.4 

i=l 

R = [2.8 ,2.3 ,2.4] 

Step 5. Ranking the preference order. 

Ranking of the alternatives would be created in the descending order of 
the valuesk (vA calculated in the fifth step. So when the fifth step in the 
calculation of the evaluation of the candidate cars (alternatives) from small to 
large k (v 2 ) < fc(v 3 ) < k(v t ), the order form is realized in the form of ranking 
v 2 < v 3 < v x . In other words, the most suitable car appears to be v x . 

3. Conclusion 

We proposed a new method for selection in this paper. After verifying the 
accumulations on different situations it can be observed that the new method is 
quite simple to use and significant for accumulation. Also, it has less no of 
calculations than the original TOPSIS. In this method we deal with the 
indeterminate and fuzzy or ambiguous values. We get same result through soft set 
on TOPSIS and with fuzzy soft set on TOPSIS. 
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Abstract 

Melia azedarach L. belongs to one of the most versatile medicinal plants family 
meliaceae (mahogany) which has great attraction for researchers. The plant was 
selected for research because it was one of the least explored members. The 
presence of saponin, alkaloids, tannins and flavonoids in the leaves extracts of 
plant indicated its medicinal value. These compounds have pharmacological 
effects against cancer, viral and malarial infections that are one of the main 
causes of deaths. With passage of time most of bacterial strains develop 
resistance against traditional medicines so they are needed to be upgraded or 
replaced. There is a need to develop antimicrobial agents with more effectiveness 
and minimum side effects. There are some reports from last two decades that 
Melia azedarach is a potential source of novel antibodies. Its extracts have both 
antioxidant and antimicrobial activities. Powdered leaves of M. azedarach were 
extracted with methanol and extract was preliminary examined by phytochemical 
tests and thin layer chromatography. The different concentrations of extracts 
showed good antibacterial activities against three pathogenic bacterial strains 
Staphylococcus aureus, Escherichia coli and Bacillus thuringiensis. The results 
indicated that M. azedarach L. could be an effective source of herbal medicines 
against infectious diseases. 

Keywords: meliaceae, melia azedarach L, antibacterial activity. 

Introduction 

The medicinal plants are major source of raw materials for treatment of various 
diseases in Ayurveda, Siddah and Unani traditional systems. In the modem 
medicine system the 25% of drugs are extracted from the plants. Azadirachta 
indica (former Melia indica) is widely used in herbal medicinal system. Melia 
azedarach L. (M. azedarach), although a close family relative, has been ignored 
in indigenous and modern system of medicine. The classification of two plants is 
presented in Table 1 show their relevance with each other. The plant is native to 
Asia and is widely spread and naturalized in most of the subtropical and tropics 
parts of the world [1-4]. Deciduous tree of M. azedarach is small to medium in 
size. Maximum height is 15m and diameter is about 60cm.The plant has 
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characteristic dark green dense and spreading crowned with dark brown bark. 
Dark green, pale leaves are short stalked, hairless and thin and white flowers 
having purple stripes and fragrance ( 1,2). It has yellow round, smooth and fleshy 
berries with four to five seeds [5]. 


Table 1. 

Classification of M. azedarach and A. indica 


Category 

Melia azedarach L. 
(dhraik) 

Azadirachta Indica( neem) 

Kingdom 

Plantea 

Plantea 

Bionomial Name 

Melia azedarach L. 

Azadirachta Indica 

Division 

Magnoliophyta 

Magnoliophyta 

Class 

Magnoliopsida 

Magnoliopsida 

Order 

Sapindales 

Sapindales 

Family 

Maliaceae 

Maliaceae 

Genus 

Melia 

Azadirachta 

Species 

M.azedarach 

M. Indica 

Synonyms 

Melia australis sweet 

Melia indica 


The extracts of leaves, fruits have antioxidant, analgesic and antimicrobial 
activities. The plants also showed antifungal, antimalarial & cytotoxic activities. 
Alkaloids, flavonoids, glycosides, saponins & tannins present in the methanolic 
extract of plant [5]. It has been used as a bitter tonic, fuel, an anthelminthic, 
astringent and an antiseptic agent [6]. It has been shown curing properties for 
diseases such as rheumatism, leprosy, rashes etc [7]. It is revealed from a 
literature survey that limonoides, limonin, nomilin and obacunone have been 
discovered [8]. The oil of M. azedarach contains sulfur containing limonoids 
which has been used in soap & cosmetic industries [9]. 

Figure 1. Leaves of Melia Azedarach Linn 



Figure2. Flowers of Melia Azedarach Linn 
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2. Materials 

2.1. Collection and Identification of Plant Material 

The plant material of M. azedarach L. was collected from garden of G.C. 
University of Lahore. Taxonomist Prof. Dr. Zaheer-Ud-Din Khan of Govt. 
College University Lahore identified the plants. Further the plant specimen was 
deposited at herbarium of G.C University. 

Family : Meliaceae 

Species : Melia Azedarach Linn. 

Voucher No.: G.C. Herb Bot. 2908 

2.2. Drying Of Plant Material 

The fresh plants leaves were placed in clean stainless tray and kept under 
shade at temperature of 20°C to 25°C for one week. After drying the leaves were 
powdered by mechanical means. 

3. Methodology 

3.1. Soxhlet Extraction 

An average of 40g (500g total) of dried and powdered plant leaves were 
taken at a time in the thimble and put into the Soxhlet apparatus. Plant material 
was extracted with 350ml of methanol for 4 to 6 hours. The crude extract obtained 
was filtered with filter paper and the filtrate was then analyzed further. The 
extract was concentrated in a rotary evaporator for complete removal of methanol 
[7]. 

3.2. Detection of Secondary Metabolites in Crude Plant Extract 

For detection of secondary metabolites, 4g of powdered plant leaves were 
extracted with methanol (50ml), boiled it for 5 minutes, cooled and filtered. In the 
filtrate added small amount of Lead acetate to make solution. Filtered and divided 
in several parts. Different tests were applied for the Tannins, Flavonoids, 
Alkaloids and Saponins according to standard procedures [10]. 

3.3. Thin Layer Chromatography 

The TLC was carried out with BAW (Butanol, Acetic Acid CH 3 COOH, 
H 2 O in ratio 4:1:5) as solvent. After drying, the chromatogram was observed 
under UV lamp. The wave length 366 nm was used for observation. Distribution 
of components was good and different colored spots were observed. Pattern of 
colored spots was recorded and Rf value of each spot was calculated. 
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3.4. Antibacterial activity 

The antibacterial investigation was carried out by well diffusion method 
against three pathogenic bacterial strains. Bacillus thuringiensis (Ref. # BL- Bt6), 
Staphylococcus aureus (ATCC # 6633) and Escherichia coli (ATCC # 25922). 
The concentrated methanolic extract of leaves of M. azedarach L was dissolved in 
distilled water and different concentrations were prepared. 

Nutrient broth was dissolved in distilled water and heated. Adjust the PH 
to 7.4 and then sterilized it in an autoclave for 25 minutes. Added some agar in it 
and heated again. Bacterial cultures were taken from stock slants. The bacterial 
cultures were incubated at 37°C for one day. These cultures were then added to 
conical flask which consists of freshly prepared nutrient broth and placed in a 
shaker at 37°C and incubated for one day. The nutrient agar in molten form was 
poured as a basal layer in petri dishes. Plates were inoculated with respective 
organism. After solidification, wells were bore which were then fill with control 
and extracts of different concentrations, and incubated again for 24 hours at 37 
°C. The diameters of zones of inhibition around the wells were recorded. The 
antibacterial activity of plant extracts were estimated by evaluation of Minimum 
inhibitory concentration (MIC) [11-13]. 

4. Results and Discussions 


4.1. Phytochemical Analysis 

Phytochemical analysis of M. azedarach methanol leaves extracts was 
carried out to investigate the absence or presence of alkaloids, flavonoids, 
saponins, glycosides and tannins in leaves and the results are presented in the 
Table2. The results supported the previous studies that plant contained Alkaloids, 
Tannins, Steroids Saponins [14]. 

Test results revealed that plant leaves contain bioactive compounds like 
Alkaloids, Steroids, Saponins, Tannins, and flavonoids. These compounds are 
responsible for the antimicrobial efficacy of plant. However the sample extract 
did not show the presence of Anthraquinones. 


Table2. 

Phytochemical analysis of methanolic extract of leaves 


Experiment Test for Alkaloids: Observation 


Mayer’s test 
Wagner’s test 
Hagner’s test 
Dragendorff test 

Test for Tannins: 


Yellow ppt. 

Reddish Brown ppt. 
Yellow ppt. 

Cloudy Appearance 


Inference _ 

Alkaloid present 
Alkaloid present 
Alkaloid present 
Alkaloid present 



SCIENTIFIC INQUIRY AND REVIEW 


Issue: I 


2017 


22 


Volume: I 







Exploration of Antibacterial Potential of Melia Azedarach L. 


5ml extract + 6 drop FeCb and 
allowed to stand for 5 minutes 
Test for Flavonoids: 

Plant extracts+2 ml AICI 3 . 

Test for Steroids: 
Leiberman-Burchard test 
Test for Saponnins: 

Shook the aqueous plant 
vigorously 

Test for Anthraquinones: 

Plant extract + few drops of N, 
V-dimethylean aline 


Bluish Green ppt. 

Yellow Fluorescence 
Violet Color 
Persistent Froth 

No red color 


Tannins present 

Flavonoids present 
Steroids present 
Saponnins present 

Anthraquinones absent 


4.2. Thin Layer Chromatography 

TLC on prepared TLC card (0.2 mm) was used for the initial investigation 
of bioactive compounds present in methanolic leaves extracts of M. azedarach. 
The chromatogram showed three differently colored spots under UV lamp 
(Table3). The Rf values of spots were compared with the standard data from 
Harborne [15] for general indication of classes of compounds present. The Rf 
values indicated the existence of Anthocyanidin-3-glycosides, Anthocyanidin 3, 
5-diglycosides or isoquercetin and rutinosides. 


Table3. 

TLC of the methanolic extract using BAW 


No. of 
spot 

Color in UV light 
(366nm) 

Rf value 

Expected compounds 

1 

Orange 

0.674 

Anthocyanidin-3-glycosides 

2 

Red 

0.584 

Anthocyanidin 3,5-diglycosides or 
isoquercetin 

3 

Green 

0.426 

Rutinosides 


4.4. Antibacterial Activity 


M. azedarach L was reported to have antibacterial properties against 
bacterial strains other than those used in this study [16,17]. The antibacterial 
activity of M. azedarach L was investigated by using three concentrations of 
concentrated methanolic leaves extract in methanol as 100 mg/ml (D2), 50 mg/ml 
(D2), and 10 mg/ml (Dl). Control used was Streptomycin for whom average 
diameter of zone of inhibition was 27.5mm. The extracts showed significant 
inhibition against Gram positive as well as Gram negative bacteria as shown in 
Table4. 
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The maximum inhibition was observed by concentration D1 in case of 
Escherichia coli while Bacillus THURINGIENSIS AND Staphylococcus aureus 
gave maximum reduction on concentration D2. This indicated the different 
response of Gram positive and Gram negative bacteria against bioactive 
compounds. 

Table4. 

Antibacterial activity of M. azedarach leaves extract 


Bacterial Strain Type of Bacteria Diameter of zone of inhibition(mm) 

Methanol Dl a D2 b D3 C 


Bacillus thuringiensis 

Gram positive 

0 

27 

28 

11 

Staphylococcus aureus 

Gram positive 

0 

27 

31 

17 

Escherichia coli 

Gram negative 

0 

39 

32 

12 


a Dl showed concentration 100 mg/ml 
b D2 showed concentration 50 mg/ml 
c D3 showed concentration 50 mg/ml 


A Comparison of MIC value of three bacteria is demonstrated in Figure3. 
It showed that methanolic extract of M. azedarach is most potent against 
Staphylococcus aureus. Therefore in view of previous reports and present study, 
M. azedarach was found to be much potent and having bioactive compounds [18- 
21 ]. 

20 



B. thuringiensis E. coli S. aureus 


Figure3. Comparison of MIC of three bacteria 

5. Conclusion 

The M. azedarach is famous for its important alkaloid content. 
Phytochemical tests and TLC showed that the plant leaves extract contains large 
variety of secondary metabolites such as alkaloids, flavonoids, saponins, tannins 
and glycosides. Due to presence of such compounds, the plant exhibited 
significant antimicrobial activity against pathogenic bacteria. The leaves extract 
of plant showed high antibacterial potential and gave significant zones of 
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inhibition. The study indicated that M. azedarach can contribute well in the herbal 
medicinal system so it can be concluded that there is a need to explore and accept 
its medicinal value. The bioactivity spectrum of plant against pathogenic bacteria 
is worth noticing. Since the bioactive compounds present in the leaves extract of 
plant are responsible for the significant antimicrobial activity against the tested 
microorganisms so this plant need to be explored further for its medical potency. 
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Abstract 

By considering density functioned theory (DFT) in terms of ab-initio investigation, 
we have explored the structured, electronic and magnetic properties of cubic 
CeCrOj for the first time. In order to determine the structured steibility of cubic 
CeCrO 3 compound, we optimized the structure of CeCrO 3 in non-magnetic (NM), 
ferromagnetic (FM) anel Anti-ferromagnetic (AFM) phases by using PBE 
generalized gradient approximation (GGA) functioned to find the exchange- 
correlation potential. From structural optimization, the FM phase of CeCrO 3 is 
observed to be stable. For computing electronic and magnetic properties, the 
lately advanced modified Becke and Johnson loced (spin) density approximation 
(mBJLDA) is used. Calculated band structures and density of states plots with an 
integer magnetic moment of 4 pB and reveal half-metallic character. In addition, 
s-d exchange constants (Noa) and p-cl exchange constant (Nofi) are determined, 
which are in agreement with a distinctive magneto-optical experiment. 

Keywords: ab-initio calculations, cubic perovskite, half-metallic ferromagnetic, 
magnetic properties 

Introduction 

Magneto-electric (ME) characteristic in the multi-ferroics, now a day’s has 
become of vigorous importance, due to the existence of correlated ordering 
parameter of electric and magnetic components that gets so much attention for 
famous existing application such as sensor of magnetic field, memory elements 
for multiple state, spintronic and multi-range microwave devices etc [1-3]. Detail 
information of crystal and magnetic structure is matter of great concern for useful 
application. Recently, the rare-earth ortho-ferrites have been reported like 
GdFeCL [5] and DyFeCL [4] have been reported for ferro-electricity as well as 
ME coupling effects. It is very difficult to perform polarization measurement 
caused by high leakage current due to high Neel temperature of RFeCL 
(TN Fe =620-740K). At the same time as perovskite chromite RC 1 O 3 exhibits 
magnetic properties at lower TN values (110-280K).RCrO3 (R=Sm- Gd- Tb- Er- 
Tm- and Y) [6,7] was reported to show signs of a fairly great electric polarization 
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(0.2-0.8 pC/cm 2 ), initially, at somewhat high temperatures equivalent to the Tn of 
the Cr sub-system. Additionally, LaCrCL, CeCrCL have the maximum Tn in 
RCrCL compounds, which is advantageous for device applications at room 
temperature [8], 

Perovskite oxides get so much attention due to its effective use in gas 
separation membranes, solid fuel cell and piezoelectric etc. [9-14]. Lanthanide 
doped perovskite type oxides, for example LnMeCL (Ln: lanthanides, Me: 
transition metals), have been accustomed for functional inorganic materials 
having a inclusive diversity of applications for alkaline fuel cells electrodes [15], 
gas ions sensors [16] and catalysts for fast and complete oxidation/reduction of 
CO, NO and other hydrocarbons [17]. Various fascinating physical properties 
such as structural, electronic, optical and magnetic properties are inter-dependent 
in transition-metal oxides [18, 19]. These materials are predictable for spintronics 
devices effectively. 

The association of magnetic and electronic played a role to develop the 
research in area of spintronics. Half-metallic ferromagnetism has a significant part 
because of its spin- polarization at the Fermi-level which is necessary for the 
better performance of spintronic applications [20, 21]. All these points force us to 
explore half-metallic ferromagnets with important magnetic moment that well- 
matched with existing semiconductor technology. Material similar to half-metal 
have a unique property to act as conductor in one direction of spin and insulator 
for opposite direction, therefore, it is very suitable for device applications and get 
attention for researcher that are working in the area of spintronic devices [22-24]. 
Some major examples are magnetic disk drives, magnetic tunnel junctions, 
magnetic hybrid technology for CMOS and magnetic sensor [25, 26], non-volatile 
magnetic random access memories (MRAM) [27, 28]. By getting research 
inspiration from de Groot et al. work, that explicate the insight of half-metallicity 
to compute band structure using the half-Heusler alloys NiMnSb [29]. Several 
research groups perform numerous experimental and computational studies on the 
HM ferromagnets, and many HM materials have been predicted and experimental 
verified [30]. 

CeCr 03 belongs to Pm-3m (No. 221) space group and have cubic crystal 
structure. The atoms arrangement are as that Cr ions are placed at center of unit 
cell and coordinated with 6 oxygen ions, Ce ions are distributed at the corner of 
the cell and oxygen at center of the faces of unit cell. We have applied mBJ 
scheme, as introduced by Becker and Johnson (BJ) [31], as it can properly find 
electronic and magnetic characteristics. The structural properties of aforesaid 
crystals are calculated at ground state and equated with the prevalent theoretical 
and experimental data. 
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2. Method of Calculations 

In this study, predicted results were obtained by carrying the Density 
Functional Theory (DFT), which is quantum mechanical approach that successful 
in predicting fundamental properties of compounds and alloys in terms of 
semiconducting trend. We employed DFT based full potential linearized 
augmented plane wave plus local orbital (FP-LAPW+lo) method within the frame 
work of Wien2K code [32]. In order to determine the ground state factors like 
lattice constant and bulk modulus, we used generalized gradient approximation 
(GGA) functional by considering the exchange-correlation potential suggested by 
Perdew, Burke, and Emzerhof (PBE) [33]. Whereas, recently developed modified 
Becke-Johnson local density approximation functional (mBJLDA) [31] were used 
to analyze the magnetic and electronic properties. The aim of using mBJLDA 
potential for electronic properties is because of that exploring the improved 
predicting bandgap as associated with standard LDA [34] or GGA [33]. 

In FP-LAPW+lo method, ion cores inside non-overlapping spheres and a 
region of constant potential (interstitial region) are considered to the region 
between the spheres. In interstitial region, a plane wave expansion is used, 
whereas basis functions, potential and charge density were prolonged as 
arrangements of spherical harmonic functions. The value of lmax=10 in muffin- 
tin spheres for charge density and non-spherical potential was accomplished. For 
energy merging, basis function expand upto RmtX K max =8 (in the plane wave 
extension Rmt represent the minimum sphere radius and Kmax the amount of the 
largest K vector). A mesh of 56 k-points was used in the irreducible part of the 
Brillouin zone (BZ) for structural, magnetic and electronic properties, which 
certifies the convergence, are 3000 k-points. In addition, the charge density was 
Fourier prolonged up to G ma x = 16. For cubic CcCrO i perovskites, the Rmt values 
were elected to be 2.5, 1.72 and 1.87a.u. (atomic units) for Ce, Cr and O 
correspondingly. For energy convergence, the calculations of self-consistent were 
performed iteratively, when the total energy of the system is steady within 0.01 
mRy. 

3. Results and Discussion 

3.1. Structural properties 

To understand DFT based cubic CeCr 03 , structure stability is check by 
performing structural optimization in NM, FM and AFM phases. By using GGA- 
PBE scheme, optimization is done by minimizing the total energy with reverence 
to unit cell volume in each phase. From computed results (see Figure 1), the total 
energy difference between these two as AE =Enm - Efm and AE =Eafm - Efm . The 
positive of AE confirmed that CeCr03 is stable in FM phase (see Table 1). In 
ABO3 oxides with B=Mn, Cr, Fe, Ni and Co, most of such oxides have FM 
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character. Similar, FM character is also proved in our study. Therefore, in first 
step, we performed optimization in FM phase to compute the lattice constant a( A) 
and bulk modulus B for CeCrCFare shown in Tablel. 

The calculated tolerance factors for CeCrCB are mentioned in Tablel. Our 
calculated value of tolerance factor is in adjacent covenant with the calculated 
results [35, 36]. In cubic perovskite, the tolerance factor lies between 0.93 and 
1.02 [36] and our measured values employed in this range, illuminating the cubic 
structure of CeCr03 compounds. The bond lengths are measured between various 
atoms of the CeCr 03 and also listed in Tablel. The tolerance factor can be 
calculated using bond lengths by using the consequent formula: 
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Figure 1. The computed total minimum energy versus unit cell volume in 
non-magnetic, ferro-magnetic and anti-ferromagnetic cubic CeCrC >3 compound 


Tablel. 

o 

Calculated lattice parameters a(A), Bulk moduli B(GPa), Tolerance factor, Bond 
length, bandgap Eg (eV), Half-metallic Ehm(cV), magnetic moments (pb) and 
exchange constant parameters of FM cubic CeCrCE 


Parameter 


CeCr0 3 


ao(A) 

3.877 

Bo(GPa) 

183.81 

Tolerance factor 

0.999 

Bond length Cr-0 

1.9373 

Bond length Ce-0 

2.7397 

Bond length Cr-Ce 

3.355 

Eg (eV) 

2.89 

Ehm(cV) 

0.38 
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Total (pb) 

4.0004 

Cr (p B ) 

2.5283 

Ce (p B ) 

0.9831 

0(p B ) 

0.0642 

Ax (d) 

4.25 

Ax (pd) 

3.14 

AEc (eV) 

0.42 

AE V (eV) 

2.70 

N 0 a 

0.33 

Nop 

2.14 


3.2. Electronic properties 

Electronic band structure with mBJ DFT studies was discussed in this 
electronic part. Energy eigen values obtained with the help of KS equation for 
Electronic band structure. In the magnetic properties of CeCrCE show match of 
band structure due to spin up (|) and down (j) orientation of FM CeCrCE 
polarized band structure in Figure2. In the analysis of spin-up state band structure, 
it will be seen that the valance band (VB) maxima and conduction band (CB) 
minima both found at point M of Briliion Zone (see in Figure2). For above 
analysis, it is seen that up spin state of CeCrCE demonstrate the FM 
semiconductors behavior. The VB maxima cross the Fermi level (EF) for spin 
down channel in CeCrCF, which describe the half-metallic manners is there for 
spin down (j) state. To disclose the starting point of density of states, partial 
density of state for CeCrCF is calculated and shown in the Figure3. 




Figure2. The Calculated Spin Polarized Ferromagnetic Band Structures for Cubic 
Cecro 3 Compound with Mbj Potential 
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The (a) for majority spin (f) and the (b) for minority spin Q) 

It understood able the Fermi level Ef is cross for spin up, total DOS, on 
the other hand the down channel form band gap at Ef, resultant the charges make 
a complete spin polarization and form this compound to utilize it for spintronics 
devices. For compound CeCr03 and by means of mBJ 2 p states and 3d states of O 
and Cr contributed primarily in the region of VB involving -4 eV to Fermi level. 
If we deeply analyzed DOS from Ef to 1 eV, we observed the 4/ orbitals of Ce 
contribution in DOS, after active over 3d states of Cr is clearly depicted. The shift 
in these states to-words high energies and the contribution of 4/ states of Ce 
increments in the region of Fermi level clues to observed p-type conductivity and 
100% spin polarization. 



Figure3. The calculated (a) total density of states (DOS), along partial 
DOS of (b) Ce, (c) Cr and (d) O by using mBJ for CeGEE compound. 

3.3. Magnetic Properties 

We obtained magnetic moment CeCr03 from mBJ and listed in the 
Table 1. Although, the calculated value is to some extent lesser than the theory 
value (po-=2.52 pB). The fact of the smaller value is that the electrons in 3d are 



SCIENTIFIC INQUIRY AND REVIEW 


Issue: I 


2017 


32 


Volume: I 















































DFT-Mbj Study of Electronic and Magnetic Properties of Cubic Cecro Compound 


not totally confined but hybridized with oxygen of its 2 p states. We have been 
also calculated magnetic moments of Ce ions. From the calculation we suggest 
that the total magnetic moment is predominantly because of Cr atoms and the Mce 
may be negligible. The half metallic characteristic of materials is verified by the 
numeral value of total magnetic moment. In our verdict the dual exchange method 
is liable for ferromagnetism seen in CeCrOs cubic perovskites and has indirect 
exchange interface among transition metals and rare-earths by the mean of anion 
O. 


In addition splitting energy A x (d) values of spin exchange that describes 
the role of 3d states of transition metal in exchange method is studied for Cr, and 
listed in the Table 1. Most exciting calculated parameters from p-d and s-d 
coupling are swap coefficients Noa, NoP, correspondingly [37], that decide the 
swap interaction among TM c/-state as well as charge carriers (holes and electrons 
in the valance band and conduction band respectively). The swap coefficients 
values for CeCrCb are shown in Table 1. Additional, the value of Noa has lower 
than NoP indicate that the s-d contact at CB minima is greatly feebler than p-d 
interaction at VB maxima, that may be evidence of ferromagnetic behavior 
present in the given compound. 

4. Conclusion 

In conclusion, we have perceived structural, electronic and magnetic 
properties of FM CeCr03 compound by means of ab-initio calculations. To verify 
the stability of CeCrC >3 compound, we have optimized the structure in PM, FM 
and AFM phases and have computed their total energy differences (AEi and AE 2 ). 
The calculated value of AEi and AE 2 are positive, which confirm that CeCrCb is 
stable in FM phase. By analyzing band structure and density of state plots, we 
find that CeCrC >3 is a half-metallic ferromagnet, while predicted value of total 
magnetic moments is 4 pB. In addition, Cr 3d (unfilled) state due to p-d 
hybridization results in a fall of magnetic moment of Cr ions and identical of 
magnetic moment in the Ce and O (nonmagnetic) ions. Furthermore, calculated 
exchange constant Noa and NoP indicate a lower value of Noa than NoP, which 
represents that spin-down state is more operative, due to p-d interaction at valance 
band maxima and ferromagnetism is confirmed by their strong hybridization. 
Furthermore, double-exchange mechanism is used to discuss the origin of the 
ferromagnetism in CeCr03. 
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Abstract 

In this communication, systematic first principles calculation has been scrutinize 
to evaluate bonding nature, structured, electronic, and optical properties of 
RbHgFs. The findings are based on total energy calculations where Khon Sham 
(KS) equation is solved by means of density functioned theory (FP-LAPW) 
method. Optimization of structured parameters is done with veiriety of 
approximations, which corroborates through comparison with available 
experimented data. Assessment of beind profile through GGA plus Trans-Blaha 
moelifieel Becke-Johnson (TB-mBJ) potential highlights underestimation of 
beinelgap with traditional Generalized Gradient approximations. Specific 
contribution of particular states on electronic properties is investigated by means 
of total anel partial density of states while contour maps of electron density are 
used to sightsee bonding character and it is evaluated that emphasized compound 
is (M-r) indirect bandgap material with mixed ionic and covalent bonding 
character. Aclditioncdly attention is paid to absorption and reflection spectra of 
RbHgFfluoroperovskite by reconnoitering optical properties, which shows 
extensive absorption and reflection in high frequency regions. Expectantly, 
current study would benchmark various quantum mechanical effects, which must 
be taken into account to understand and utilize RbHgFj in fabricating practical 
devices. 

Keywords: first- principles study, fluorine based perovskites, electronic property, 
optical property 

Introduction 

The knowledge of physical properties of materials is always a principal field of 
curiosity but in today’s hi-tech era trust of potential investors can be strengthened 
by engrossment of researchers in competition towards higher proficiencies. To 
fulfill this need material scientists & researchers struggling hard for economical 
and resourceful materials. The class of compounds having ABL3 stoichiometry is 
known as fluoroperovskites. Here A is usually alkali, alkaline or rare earth metals 
while B is supposed to be transition, post transition and non-transition metals 
although anion is represented by X that are oxides and halides [1], 
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In this study emphasize is given to Rubidium based mercury 
fluoroperovskite which own great technological importance in fabricating low 
birefringent excellent lens materials, transparent coating devices, opto-electronic 
& optical pathways, photovoltaic applications, ionic conducting, luminescence 
capacitor, UV detectors as well as Light Emitting Diodes (LED). Particularly 
RbHgF 3 is a superior choice for the high-class lens materials because it does not 
suffer with birefringence, which can make design of lenses problematic [2-10], 
Comprehensive experimental studies are available on their interesting structural 
properties which is being confirmed that under consideration compound 
crystallizes in cubic structure and it do not unveil any phase transition under 
varies pressure and temperature [11-12], 

From above literature this idea can be established that focused material is 
technically sound which motivate us to inquire it in detail. Our main goal is to 
improve trends of all physical properties concurrently. This paper is schematized 
in the following sections where Section 1, have already described some 
introductory details about materials and their applications. Section 2, is devoted to 
method and computational detail of calculations. Section 3, present results, and 
discussions of structural and opto-electronic properties. Finally, all results are 
compared with previous studies where data is available. At the end, we enlighten 
future prospective of this research. 

2. Computational Details 

The present first principles study have been carried out on the basis of 
density functional theory (DFT) which is implemented in wien2k code to solve 
Kohn Sham equation within Full Potential L inearized Augmented Plane Wave 
(FP-LAPW) method. This method is one of the best methods for appropriate 
computation of electronic states and optical response of various crystalline solids. 
For structural optimization the exchange correlation approximation is treated 
with Perdew-Burke-Emzerhof Generalized Gradient Approximation (PBE-GGA), 
Wu Cohen Generalized Gradient Approximation (WC-GGA), PBEsol-GGA and 
Local Density Approximation (LDA) while for opto-electronic properties recently 
bugged GGA plus Trans-Blaha modified Becke-Johnson (TB-mBJ) potential is 
employed [13-16]. 

In order to enhance accuracy in calculations valence electrons are treated 
semi-relativistically and core electrons are treated fully relativistically. 
Convergence in basis size is achieved with a cut-off RMT.Kmax = 8.0, which is 
product of smallest muffin-tin radius RMT times the largest plane wave vector 
Kmax. Brillouin zone (BZ) integration is done with 56 K point using modified 
form of tetrahedron method [17]. The calculations are self-consistently converged 
when total energy and charge is stable within 0.00ImRy and 0.01 m respectively. 
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Furthermore, to calculate validate optical properties denser mesh of Brillouin 
zone (BZ) sampling is done with enormous number of K points. 

Table 1. 

Comparison of Present calculation with previous experimental and theoretical 
values for lattice constants (a 0 ), ground state energies (E 0 ), bulk modulus (B 0 ) and 
its pressure derivative (B p ) of RbHgF 3 compound. 


Compound 

RbHgF 3 

Present work 

Present work 

Present work 

Present 

work 

LDA 

Experi 

mental 

work 

Other 
theoretic 
al work 

PBE-GGA 

WC-GGA 

PBEsol- 

GGA 

a„ (A) 

4.60 

4.57 

4.53 

4.49 

4.47 a 

4.46 b 

E 0 (Ry) 

-45854.51 

-45854.43 

-45854.40 

-45854.39 



B„ (GPa) 

48.84 

49.01 

49.39 

49.81 

48.32 a 


Bp(GPa) 

5.61 

5.58 

5.53 

5.51 




a) Ref. [22] (Experimental work) b) Ref.[23] (Other theoretical work) 

3. Results and Discussion 

3.1. Structural Properties 

Ternary fluoroperovskite RbHgF 3 crystallizes itself in cubic type of 
structure having space group Pm3m (#221) while its unit cell occupy one 
molecule. Figure 1 illustrates the sites of Wyckoff coordinates which are situated 
at la (0 0 0), lb (1/2,1/2,0), 3c (0, 1/2, 1/2) for Rb, Hg and F respectively. This 
subsection is dedicated to calculate structural properties of RbHgF3 via energy 
minimization process in which total energy varies as a function of equilibrium cell 
volume (Vo). Accurate total energy versus volume curve is fitted with equation of 
state (EOS) developed by Murnaghan [18] as shown in Figure2. 

Equilibrium lattice constant (a 0 ), ground stata energy (E 0 ), bulk modulus 
(B 0 ), and pressure derivative of bulk modulus (Bp) are calculated with four 
different exchange and correlation schemes namely PBE-GGA, WC-GGA, LDA 
and PBEsol-GGA. However from Table 1 it can be observed that lattice constant 
computed by DFT slightly overvalue as compared to experimental data that can 
be associated due to use of traditional DFT schemes. The value of bulk modulus 
represents good crystal rigidity in RbHgFsfluoroperovskite compound. 
Furthermore, Table 1 depicts contrary relation between lattice constant and bulk 
modulus in accordance with the trend of other fluoroperovskite compounds [19- 
21]. As a result all structural parameters are in reasonable agreement with 
previous theoretical and existing experimental data. 
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Figurel. Crystal Figure2. Variation of total 

structures of RbHgF 3 energy as a function of unit 

cell volume for RbHgF 3 

3.2. Electronic Properties 


The electronic properties of RbHgF 3 fluoroperovskite compound is 
described in terms of calculating energy band structure and total as well as partial 
density of states (TDOS and PDOS) whereas bonding nature is evaluated with aid 
of electron density plots. In this subsection band gap is computed with GGA plus 
Trans-Blaha modified Becke-Johnson (TB-mBJ) potential in order to avoid 
underestimation of well-known exchange and correlation scheme of LDA and 
GGA [24-27], To clarify this concept bandstructure comparison is made with 
PBE-GGA scheme to prove this underestimation. The resultant energy 
bandstructure in high symmetry directions are shown in Figure3 with PBE-GGA 
and mBJ potential. It can be analyzed from s that overall trend of band dispersion 
curves are almost same and conduction band minimum (CBM) lies at T symmetry 
point of brillouin zone (BZ) whereas valence band maximum (VBM) is located at 
M point brillouin zone (BZ) which reveals (M-T) indirect bandgap of 3.2 eV and 
1.0 eV from mBJ and PBE-GGA schemes respectively. Unfortunately, there is 
lack of experimental bandgap data to make a reasonable comparison. However, 
the compound can work well in ultraviolet region of electromagnetic spectrum 
because the materials with band gaps larger than 3.1 eV work well in the 
ultraviolet region of the spectrum [28]. 

Energy density distribution of varied states is observed in terms of total 
and partial density of states (TDOS & PDOS). According to Figure4-5 it is 
obvious that TDOS and PDOS of RbHgF 3 fluoroperovskite compound can be 
splitted into varied energy regions ranging from -10 eV to 15 eV. A narrow sharp 
peak is observed at -10 eV due to Rb-4p states. The upper region of valence band 
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from Fermi-level to -6.1 eV is due to overlapping of Hg-3d and F-2p states. 
While above Fermi level, lower part of conduction band is occupied by Hg-4s and 
upper part of conduction band is filled by Rb-4d states. For the sake of precision 
in density of states, we omit some distorted peaks which can hinder justification 
of electronic states. 




Figured. Comparison of calculated band structure in high symmetry directions with PBE- 

GGA & mBj potential 
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Figure4-5 Calculated total and partial density of states (TDOS & PDOS) 
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(I) (100) plane 


(II) (110) plane 


Figure6. Calculated mBJ total two and three dimensional electronic charge 
densities for (I) in the (100) plane and (II) in the (110) plane. 

The maps of contour plots in terms of electron density communicate a 
crucial role in explaining nature of chemical bond in crystalline materials [29], 
The dispersion curves of contour maps are calculated along (100) and (110) 
planes in 2D alongwith 3D as quoted in Figure6. The strong ionic nature of Rb-F 
bond is evident from perfect spherical charge distribution between cation of 
rubidium and anion of fluorine. The charge is transferred between Rbcation and F 
anion due to large electronegativity difference in them. On the other hand, 
uniform distribution is observed between Hg cation and F anion which divulges 
covalent nature in HgFaoctahedra. In fact above analysis synchronizes well with 
results of DOS Figure4-5 where p-d hybridization is maximum between Hg-3d 
and F-2p states. Conclusively our results authorize mixed ionic as well as 
covalent bonding nature in RbHgF 3 . Furthermore in Refs. [20,21] similar results 
of bonding type is observed for other fluoroperovskite compounds. 

In order to expose internal behavior of emphasized compound the tool of 
optical analysis is employed. Fundamental optical responses include imaginary 
part of dielectric function s 2 (co), real part dielectric function 
^(co). Further analysis is extended to plot optical spectra of absorption coefficient 
and reflectivity which are crucial in concern with device applications. 
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The complex dielectric function £( 0 ) can break into real and imaginary 
part according to following relation [30]: 

e(<u) = £1 (<u) + ie 2 (<u) (1) 

In equation (1) £ 2 (®) denotes imaginary part of dielectric function while 
£ 1 ( 0 ) represent real part of dielectric function. The analysis of £ 2 ( 0 ) as shown in 
Figure7 out complete response of material due to applied electromagnetic 
radiation. In this study, we highlight results of direct interband transition although 
summarizing all possible transitions originating from occupied valence band to 
unoccupied conduction band taking into account appropriate element of transition 
dipole matrix [31].The widespread peaks of £ 2 ( 0 ) follow pattern of DOS and 
band structure of the investigated compound. The threshold energy point occurs at 
5.2 eV approximately while the major peak are located at 5.9 eV which 
corresponds to transition of occupied valence band states to unoccupied 
conduction band states. These peaks are majorly due to Hg-4s states. After that 
diverse peaks are observed till 15eV which occurs due to hybridized states of Rb- 
4d alongwith some p states of Hg as well as F. 

Absorptive behavior of material is analyzed with the help of real part of 
dielectric function fjX©) as shown in Figure8. At zero frequency limit static part 
of dielectric function £ x (0) is observed at 1.88 eV. The curves of £ a fco) starts 
increasing and attains a maximum value at around 5 eV and an overall narrow 
bandgap semiconductive nature is observed. 

The plot of absorption coefficient as a function of energy depicts that 
RbHgF3 starts absorbing electromagnetic radiation at about 5.25 eV as presented 
in Figure9. This particular energy is also known as threshold point. This threshold 
point is exactly in accordance with trend of bandgap. Hereafter oscillations are 
observed with increasing absorption pattern in spectra. Investigated material starts 
absorbing effectively within 21-25 eV range while highest prominent peak is 
observed at around 21.5 eV. After the incent of optimum absorption peak it again 
going to decrease suffering trivial variations. Analysis of absorption spectra 
concludes clearly about application of this fluoroperovskite for absorption 
purposes in wide range of Ultra-Violet region of electromagnetic spectrum 
typically at about 21.5 eV. 

The calculated spectrum of reflectivity as a function of electromagnetic 
energy is shown in FigurelO. The phenomenon of reflectivity stays below 9-10% 
upto 20eV. Currently focused material starts reflecting highly and attains 
maximum value almost at 22eV that helps to understand that by following 
different regions of electromagnetic spectrum (EM) RbHgFi remains highly 
transparent in infrared (IR) and visible (V) regions which recommends that under 
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consideration compound is an upbeat candidate for efficient lenses and 
transparent coating devices. 



Figure7. Calculated imaginary part £ 2 ( 6J ) of the di-electric function for RbhgF 3 

compound 




Energy (eV) 

Figure8. Calculated real part ^(co) of the dielectric 
function for RbHgF 3 compound. 
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Energy (eV) 

Figure9. Reflectivity R(co) as a function of energy for 
RbHaF^ conmound 



Figure 10. Absorption coefficient a(co) as a function of energy for RbHgF 3 

compound 

4. Conclusion 

In this work, all electron self-consistent Full Potential- L inearized 
Augmented Plane Wave (FP-LAPW) method is used to explore structural, 
electronic, and optical properties of RbHgF 3 fluoro-perovskite compound. In order 
to inspect optimized structural parameters exchange-correlation contribution is 
treated with four different approximations. To elude underestimation of band 
profile by other DFT schemes we highlight GGA plus Trans-Blaha modified 
Becke-Johnson (TB-mBJ) potential in lieu of attaining opto-electronic results 
near to expected experimental findings. Energy band profile confirms that 
investigated material is a narrow and indirect energy bandgap (M-T) 
semicondutor. The total and partial density of state curves is used to define 
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contribution of different bands. In addition to it contour maps of electron density 
authenticates mixed covalent-ionic behavior. These results are in favorable 
agreement with previous theoretical and existing experimental data. The optical 
properties are discussed in terms of advantageous functional parameters and 
analysis is being done by interband contribution that shows RbHgFs possess wide 
range of absorption and reflection in high frequency regions and these 
characteristics make them useful for transparent coatings and fabricating low 
birefringent high-quality lenses. Conclusively limelight fluoroperovskitecan 
efficiently be used for fabricating high-class low birefringentlense material and 
photovoltaic applications. 
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Abstract 

Medicinal plants are conventionally used for the treatment of various diseases 
due to their world-wide occurrence and least side effects. Melia azedarach L. 
belongs to the family Meliaceae, is a highly significant medicinal plant. Extracts 
of M. azedarach obtained from its different parts such as seed, fruit, flower, leaf, 
and young branches are reported to exhibit antifungal, antihelmintic, nematicidal, 
diuretic, cytotoxic, antiproliferative, insecticidal and antioxidant activities. Thus 
the aim of this study was to explore the chemical profile of non-polar extract of 
M. azedarach leaves through GC-MS analysis. The identification of 
phytochemical compounds is based on molecular ion peak, base peak, and 
fragmentation pattern. GC-MS analysis of hexane extract of M. azedarach 
showed a highly complex profile, containing ketones, ethers, fatty acid 
derivatives, methyl esters, 1,3-dipalmitate, 7,8-dihydrocarpesterol, and 2- 
Undecanol. This study will be useful to explore the active components of 
medicinal plants and can validate their medicinal value. 

Keywords: Melia azedarach L., Meliaceae, GC-MS method, Phytochemical 
compound. 

Introduction 

There is an increasing demand for the herbal medicine to cure variety of diseases 
as these medicines are efficient and safe without any side effects as compared to 
synthetic drugs [1]. For this purpose, scientist, botanists, chemists, and 
pharmacists all over the world are working on herbal medicines [2]. Melia 
azedarach Linn. (M. azedarach) commonly known as chinaberry or Persian lilac 
tree is a plant species of the family Meliaceae that contain 45 genus and 750 
species. It is deciduous tree of medium sized that grows to a height of five to 
fifteen meter tall and thirty to sixty cm in diameter. It can be grown successfully 
in a wide variety of situations even in alkaline soil where other trees might fail to 
grow. It is native to Indochina, Pakistan, India, Australia, and Southeast Asia. It is 
cultivated in most of the countries located in tropics and subtropical region. This 
plant has long been known as an insecticidal and medicinal plant all over the 
world due to its world-wide availability and fewer side effects. Extracts of M. 
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azedarach obtained from its different parts such as seed, fruit, flower, leaf, and 
young branches have been used for the treatment of diabetes, malaria, intestinal 
worms, cough, nausea, vomiting and paroxysmal fever, and skin disease [3-5], 

GC-MS has been known as a powerful technique for providing metabolic 
profiling of plants [6-9], Hexane extracts of many species of Meliaceae family are 
analyzed by this technique but there has been no report on GC-MS of hexane 
extract of M. azedarach or any species of genus Melia. Non polar profile of M. 
azedarach has not been explored yet. Therefore, the aim of the study was to 
explore the chemical profile of non-polar extract of M. azedarach leaves through 
GC-MS analysis. 

2. Material and Methods 

2.1. Preparation of plant extract: The fresh leaves of plant M. azedarach was 
collected from the Botanical Garden of Government College University, Lahore 
and washed individually to remove impurities and dried under shade. Dried leaves 
were crushed into fine powder using a grinder. 

The dried plant powder were weighed and dipped in a hexane and left for seven 
days. The hexane extract was then filtered. The hexane extract was concentrated 
under reduced pressure and low temperature in a Rotary Evaporator. The semi 
solid extract were obtained. 

2.2. Phytochemical analysis of plant extract: Phytochemical analysis were 
carried out by employing standard procedures to sort out the presence of 
flavonoids, alkaloids, tannins, saponins and steroids in hexane extract of Melia 
azedarach [10, 11]. 

2.3. Preparation of sample for GC/MS analysis: A measured amount of semi 
solid was re-dissolved in //-hexane of GC grade and micro-filtered to prepare the 
samples of 5.0 mg /10 mL. 

2.4. Gas Chromatography Mass Spectrometry analysis: A GC-MS analysis 
was carried out on a Shimadzu GCMS-QP2010A system in El mode (70 eV) with 
DB-5MS capillary column (30 m x 0.25 mm i.d., film thickness: 0.25 pm, J and 
W scientific, Fulsom, CA, USA.l pL of samples were injected at 250 °C with a 
split ratio of 50/50 under electronic pressure to maintain a constant flow (0.67 
mL/min) of the helium carrier gas. The oven temperature was programmed from 
150 °C for 4 min and heated to 300 °C at a rate of 3 °C/min and kept constant at 
this temperature for 2 min. The mass spectrometer was set to scan the mass range 
40-600 amu with ion source temperature 200 °C and interface temperature was 
250 °C. Analyses were performed in triplicate with a blank run after every 
analysis. The resulting data was interpreted by using Shimadzu Lab Solution, GC- 
MS Postrun analysis software. Compounds were identified by comparing the ions 
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fragmentation pattern of mass spectra with those of known compounds stored in 
NIST 147 and NIST 127 libraries. 

3. Results and Discussion 

The genus Melia is famous for the presence of liminoids and terpenoids 
that are water-insoluble plant components [12]. Previous researches have proved 
that these compounds have diverse effect against pathogens. Moreover, non polar 
profile of M. azedarach has not been explored yet. The waxy coating of leaves 
can only be dissolved by non polar solvent because waxes are non polar in nature. 
Therefore, the leaves of M. azedarach were extracted with n-hexane in order to 
isolate the wax. 

3.1. Phytochemical analysis of plant extract: The phytochemical analysis was 
carried out and it was seen that alkaloids and flavonoids are absent in hexane 
extract of M. azedarach, whereas it contain tannins, terpenes, quinines, saponins 
and steroids as given in Table 1. 


Table 1. 

Phytochemical Analysis 


Experiment 

Observation 

Inference 

Test for Alkaloids 

Mayer’s test 

No yellow ppt. 

Alkaloids absent 

Wagner’s test 

No reddish brown ppt. 

No cloudy appearance 

Alkaloids absent 

Dragendorff test 

Bluish green ppt. 

Alkaloids absent 

Test for Tannins 

5ml extract + 6 drops FeCb 


Tannins present 

Test for Coumarins 

Yellow floresence 


plant extract+ IN NaOH+ 

Pass through Uvlight 


Coumarins present 

Test for Flavonoids 

yellow florescence is not produced 


Plant extract +2ml AlCb 


Persistent froth 

Flavonoids absent 

Test for Saponins 

Shook the aqueous plant 

vigorously 

Violet color 

Saponins present 

Test for Steroids 

Leibermann-Burchard test 

Red color 


Test for Anthraquinon 


Steroids present 

Plant extract + few drops N.N- 

dimethylanaline 


Anthraquinonses 

present 


3.2. Gas Chromatography Mass Spectrometry analysis: Gas Chromatography 
and Mass spectroscopy analysis of compounds was carried out in n-hexane leaf 
extract of M. azedarach (Figure 1). 
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Figure 1. GC chromatogram of n- hexane extract of M. azedarach leaves 

Different compounds were identified in the hexane extract of M. 
azedarach on the basis of retention time, molecular ion peak, base peak, and mass 
fragmentation pattern. Eight peaks from non-polar plant extract were identified 
as: 2-Undecanol, Methyl 4, 6-decadienyl ether, 13-Docosenoic acid, 7, 8- 
Dihydrocarpesterol, Glutaric acid, dimethyl ester, Nonanoic acid, 1, 2, 3- 
propanetriyl ester, Glycerol 2-acetate 1, 3-dipalmitate and Docosenoic acid, 1 
methyl-butyl ester as shown in Table2. 

Table2. 

Compounds from hexane extract of M. azedarach 


Name of compound 

Retention Time (mins) 

Base Peak (m/z) 

2-Undecanol 

22.39 

45 

Methyl 4, 6-decadienyl ether 

22.45 

45 

13-Docosenoic acid 

33.88 

55 

7, 8-Dihydrocarpesterol 

45.48 

414 

Glutaric acid, dimethyl ester 

47.24 

133 

Nonanoic acid, 1, 2, 3-propanetriyl ester 

52.09 

151 

Glycerol 2-acetate 1, 3-dipalmitate 

52.18 

40 

Docosenoic acid, 1 methyl-butyl ester 

54.90 

44 


2-Undecanol was eluted at 22.39 min. Base peak of the compound was 
observed at m/z 45. The fragment ions observed m/z 40, 60, 83, 97, 150. The 
fragment ion observed at m/z 83 was due to the loss of -CH2 group from 
molecular ion (M+). Methyl 4, 6-decadienyl ether was eluted at 22.45 min. Base 
peak of the compound was observed at m/z 45. The fragment ions observed m/z 
41, 56, 70, 84, 150. The fragment ion observed at m/z 70 was due to the loss of - 
CH2 group from M+. 13-Docosenoic acid was eluted at 33.88 min. Base peak of 
the compound was observed at m/z 55. The fragment ions observed m/z 69, 83, 
97, 185, 207, 320. The fragment ion observed at m/z 83 was due to the loss of - 
CH2 group from M+. 7, 8-Dihydrocarpesterol was eluted at 45.48 min. Base peak 
of the compound was observed at m/z 441. The fragment ions observed at m/z 91, 
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133, 147, 191, 208, 335, 441. The fragment ion observed at m/z 191 was due to 
the loss of -OH group from M+. Glutaric acid, dimenthyl ester was eluted at 47.24 
min. Base peak of the compound was observed at m/z 133. The fragment ions 
observed m/z 40, 69, 83, 193, 209, 281.The fragment ion observed at m/z 40 was 
due to the loss of -C2H5 group from M+. Nonanoic acid, 1, 2, 3-propanetriyl ester 
was eluted at 52.09 min. Base peak of the compound was observed at m/z 151. 
The fragment ions observed m/z 44, 68, 85, 100, 169, 207, 281, 335. The 
fragment ion observed at m/z 68 was due to the loss of -OH group from M+. 
Glycerol 2-acetate 1, 3-dipalmitate was eluted at 52.18 min. Base peak of the 
compound was observed at m/z 40. The fragment ions observed m/z 68, 100, 169, 
209, 267, 335, 441. Docosenoic acid, 1 methyl-butyl ester was eluted at 54.90 
min. Base peak of the compound was observed at m/z 44. The fragment ions 
observed m/z 77, 101, 137, 150, 208, 341. The fragment ion observed at m/z 177 
was due to the loss of -CH3 group from M+. 

Conclusion 

Traditionally, medicinal plants are largely used to cure many diseases 
owing to their least side effects. M. azedarach L. belongs to potent medicinal 
plant family Meliaceae. This was the first report on GC-MS of Hexane extract of 
M. azedarach. The results showed that non polar extract of the plant contained a 
variety of compounds such as ketones, ethers, fatty acid derivatives, methyl 
esters, 1,3-dipalmitate, 7,8-dihydrocarpesterol, and 2-Undecanol. This will 
contribute in completing the chemical profile of M. azedarach. 
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